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I.  INTRODUCTION 


It  has  been  shown  [1]  that  the  performance  of  an  adaptive  antenna 
array  Is  strongly  affected  by  the  electromagnetic  characteristics  of  the 
antenna  array.  An  Important  electromagnetic  characteristic  of  an 
antenna  array  Is  the  mutual  coupling  between  Its  elements.  In  the  above 
work,  mutual  coupling  between  the  antenna  elements  was,  however. 

Ignored,  l.e.,  the  antenna  elements  were  assumed  to  be  Isolated  from 
each  other.  In  practice,  elements  of  an  antenna  array  have  mutual 
coupling,  which  In  turn  affect  the  gain,  beamwldth,  etc.,  of  the  array. 
Mutual  coupling  becomes  particularly  significant  as  the  Interelement 
spacing  Is  decreased. 

In  this  report,  the  effect  of  mutual  coupling  on  the  performance  of 
adaptive  arrays  Is  studied.  It  Is  shown  that  the  mutual  coupling  does 
affect  the  performance  of  adaptive  arrays  and  these  effects  are 
significant  even  for  large  Interelement  spaclngs,  l.e.,  for  spacing  of 
more  than  half  a  wavelength.  The  effect  Is  rather  drastic  as  the 
Interelement  spacing  drops  below  a  half  wavelength.  In  fact,  for  a 
fixed  aperture  with  half  wavelength  spaced  elements,  the  Introduction  of 
additional  elements  can  degrade  the  array  performance.  The  failure  to 
recognize  the  presence  of  mutual  coupling  will  degrade  the  performance 
of  Applebaum  type  adaptive  arrays  more  than  that  of  LMS  arrays  since  the 
optimum  excitation  has  to  be  modified  both  In  phase  and  amplitude  to 
Include  the  changes  In  the  desired  signal  vector  due  to  the  presence  of 
mutual  coupling. 


In  section  II,  an  analytic  expression  for  the  steady  state  output 
slgnal-to-lnterference-plus-nolse  ratio  (SINR)  of  an  adaptive  array  Is 
derived.  The  expression  takes  Into  account  the  mutual  coupling  between 
the  array  elements  and  Involves  the  normalized  impedance  matrix  of  the 
array  elements.  The  expression  Is  used  to  study  the  effect  of  mutual 
coupling  on  the  performance  of  adaptive  arrays  and  It  Is  shown  that  the 
output  SINR  of  the  array  depends  upon  the  mutual  coupling  between  its 
elements.  In  section  III,  the  effect  of  mutual  coupling  on  the 
transient  performance  of  adaptive  arrays  is  studied.  It  Is  shown  that 
the  presence  of  mutual  coupling  between  the  array  elements  reduces  the 
speed  of  response  of  an  adaptive  array.  In  section  IV,  the  optimum 
excitation  to  maximize  the  output  SINR  of  Applebaum  type  adaptive  arrays 
In  the  presence  of  mutual  coupling  Is  found.  Section  V  contains  our 
conclusions. 

II.  STEADY  STATE  PERFORMANCE  OF  AN  ADAPTIVE  ARRAY  IN  THE  PRESENCE  OF 

MUTUAL  COUPLING 

The  output  SINR  of  an  adaptive  array  is  the  most  commonly  accepted 
measure  of  Its  steady  state  performance,  and  will  accordingly  be  derived 
first.  The  expression  takes  Into  account  the  mutual  coupling  between 
the  array  elements.  The  expression  will  be  used  to  study  the  effect  of 
mutual  coupling  on  the  performance  of  various  adaptive  arrays. 

The  basic  diagram  of  an  adaptive  array  is  shown  in  Figure  1.  The 
output  signal  from  each  antenna  element  Is  multiplied  by  a  complex 


Figure  1.  Basic  adaptive  array 
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weight  and  then  these  signals  are  summed  to  produce  the  array  output 
S(t).  The  weights  are  automatically  adjusted  to  optimize  the  output 
SINR  in  accordance  with  a  selected  algorithm.  To  find  an  expression  for 
the  output  SINR,  one  should  know  the  element  output  voltages.  We  will, 
therefore,  first  develop  an  expression  for  the  element  output  voltages 
when  the  mutual  coupling  Is  taken  Into  account.  These  voltages  will  be 
used  as  the  Input  signals  to  the  adaptive  processor.  The  required 
expression  can  be  obtained  by  considering  the  N-element  array  as  an  N+l 
terminal  linear,  bilateral  network  responding  to  an  outside  source  as 
shown  In  Figure  2. 

Referring  to  Figure  2,  each  port  of  the  N-element  array  Is  shown 
terminated  in  a  known  load  Impedance,  Zl.  The  array  has  as  Its  driving 
source  a  generator  with  open  circuit  voltage  Vg  and  internal  impedance 
Zg.  Using  standard  notation,  one  can  write  the  Klrchoff  relation  for 
the  N+l  terminal  network  as 


v1  =  1j  Zn  . . .  1j  Zij  . . .  1N  Zin  +  is  zls 


vi  =  1\  Zji  . . .  1j  Zjj  + . .  1  m  Zjn  +  is  Zjs  /  (2*1) 


ft 


vN  =  1 l  Zni  + 


+  1j  zNj  + 


+  In  zNN  +  is  zNs 


Further,  making  use  of  the  relationship  between  terminal  current 
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(2.2) 


If  all  the  elements  in  the  array  are  in  an  open  circuit  condition 

then 

ij  =  0  j  =  1,  2,  ....  N  , 

and  from  Equation  (2.1) 

vj  =  voj  3  zjsU  •  (2-3) 

Substituting  Equations  (2.2)  and  (2.3)  into  Equation  (2.1)  one 
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Or,  more  compactly 


Z0V  =  V0  .  (2.5) 

In  Equation  (2.5)  Z0  is  the  normalized  Impedance  matrix  and  V0 
represents  the  open  circuit  voltages  at  the  antenna  terminals.  Since  Z0 
is  nonsingular,  one  can  find  the  element  output  voltages  from  the  open 
circuit  voltages.  The  element  output  voltages  will  be  given 


It  should  be  noted  that  the  matrix  Z0  is  a  normalized  impedance 
matrix,  normalized  to  the  load  impedance.  It  acts  like  a  transformation 
matrix,  transforming  the  open  circuit  element  voltages  to  the  terminal 
voltages.  What  is  normally  assumed  in  analyzing  adaptive  antenna 
systems  is  that  the  element  spacing  Is  large  enough  so  that  the  mutual 
coupling  between  the  elements  is  small  and  consequently  the  matrix  Z0 
becomes  diagonal.  If  one  further  assumes  that  the  self -impedances  (Zii, 

1=1,2 . N)  are  equal,  the  input  signal  vector  will  be  just  the  open 

circuit  voltage  vector  multiplied  by  a  trivial  scaling  factor  Involving 
the  self  and  load  impedance  terms.  Thus  the  array  performance  will  be 
the  same  as  calculated  using  the  open  circuit  voltages  as  the  Input 
signal  to  an  adaptive  processor. 

Let  m+1  CW  signals  (one  desired  and  m  jammers)  of  the  same 
frequency  be  incident  on  the  array.  Then  the  open  circuit  voltages  at 
antenna  terminals  are  given  by 


v0  “  xd  +  l  Xik  (2.7) 

k*l 

where 

Xd  *  Ad  ej(a,ot+  ^  Ud  (2.8) 

V  *,k  *JKt  +  Mk)  »,k  • 

In  Equations  (2.8)  and  (2.9),  Ad2  i§  the  average  power  In  the 
desired  signal,  A-jk2  the  average  power  In  the  kt(l  jammer,  tug  Is  the 

carrier  frequency,  i(y  Is  the  carrier  phase  of  the  desired  signal  at  the 
coordinate  origin,  Is  the  carrier  phase  of  the  kth  jammer  at  the 
coordinate  origin  and  Ud,  U^k  are,  respectively,  the  desired  signal 
vector  and  the  k^  jammer  vector  defined  as  follows: 


fl(0d,  $d,  pd)  e^Pdl 

VV  p„)  ^2 
W  V  eJPdN 


(2.10) 


where  (6d,  <^)  defines  the  desired  signal  direction,  pd  is  the 
polarization  of  the  desired  signal,  fj(9,<l>,p)  is  the  pattern  response  of 
the  j^11  element  to  a  signal  incident  from  direction  (0,$)  with 


i 

( 


polarization  p  and  Pdj  Is  the  desired  signal  phase  at  the  element, 
measured  with  respect  to  the  coordinate  origin. 


(2.11) 

where  the  notation  Is  analogous  to  that  for  the  desired  signal  vector. 

Using  Equations  (2.6)  and  (2.7),  the  Input  signal  to  the  adaptive 
processor  will  be 

-1  m 

V  =  Z0  (Xd  ♦  I  Xik)  .  (2.12) 

k=l 

If  thermal  noise  Is  also  added  to  each  element  of  the  array  then 
the  total  input  signal  to  the  processor  will  be 

X  =  V  +  Xn 

-1  m 

=  Z0  +  l  Xikl  +  Xn  (2.13) 

k=l 

where  Xn  is  the  noise  vector  defined  as 


W  pik> eJp,kl 

W  Pik) eJp,k2 

fN ^ 9ik *  ^ik*  pik ^  ejPil<N 


Xn  =  (ni(t),  n2(t) ,  ..., 


nN(t))T 


(2.1 


In  Equation  (2.14),  T  denotes  transpose.  In  the  case  of  an 
adaptive  array,  the  signal  xj(t)  from  the  j**1  element  Is  multiplied  by  a 
complex  weight  wj(t).  The  signals  are  then  summed  to  produce  the  array 
output.  Using  the  IMS  algorithm  [2],  the  steady  state  weight  vector,  W, 
of  the  array  Is  given  by 

W  »  *-lS  (2.15) 

where  *  Is  the  covariance  matrix 

*  =  E{X*XT)  (2.16) 

and  S  Is  the  reference  correlation  vector 

S  *  E{X*R(t)}  .  (2.17) 

In  Equations  (2.16)  and  (2.17),  R(t)  Is  the  complex  reference 
signal  In  the  adaptive  array  [2,3],  *  denotes  complex  conjugate  and  E{ - } 
denotes  expectation. 

From  Equations  (2.13)  and  (2.16) 

_1  m  _1  m 

♦  -  E  {[Z„  (Xd  *  I  Xlk)  +  X„]*  [Z0  (Xd  ♦  l  X1k)  *  X„]T  } 
k*l  k*l 


_1  m  m 

(Z0  )*E  (C(Xd  +  l  Xik)  +  Z0Xn]*  [Xd  +  l  Xik)  ♦  Z0Xn]T}(Z0-l)T  . 
k=l  k=l 


(2.18) 
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Assuming  that  the  thermal  noise  voltages  from  the  array  elements 
are  gausslan  with  zero  mean  and  are  uncorrelated  with  each  other,  and 
the  carrier  phases  of  the  narrowband  signals  are  uniformly  distributed 
on  (0,  2w)  and  are  statistically  Independent  of  each  other  and  of  the 
thermal  noise  voltages,  the  covariance  matrix  *  Is  given  by 


-1  *T  m2*T  2*T  -1 

*-(Z0  )*[<?l0Z0  ♦  l  Aik  Uik  Uik  ♦  Ad  Ud  UdJ(Z0  )T 

k=l 

where  i$  the  thermal  noise  power.  From  Equation  (2.19) 


(2.19) 


♦  =  (Z0~V  o2[Z*  zj  +  "  dk  U*k  ujk  ♦  5d  Urf  ujj  (Z0“1)T 


(2.20) 


where  di  Is  the  ratio  of  the  desired  signal  power  to  the  thermal  noise 
power  and  dk  is  the  ratio  of  the  kth  jammer  power  to  the  thermal  noise 
power.  Let 

Rn  =  Z0  ZQT  +  l  dk  uik  ^ik 

k  =1  (2.21) 


♦  *  (ZqV  o2  (Rn  +  Cd  uJ)(Z0’V 


(2.22) 


Note  that  Rn  is  the  normalized  (with  respect  to  the  thermal  noise 
power)  covariance  matrix  of  the  undesired  signals  (jammer  and  the 
thermal  noise).  To  find  the  steady  state  weights  (Equation  (2.15)), 
must  be  computed.  The  following  matrix  inversion  Lemma  [41,  is  used  to 
compute 


(A  -  ozV)-1  -  A->  -  BA-'zYa-1 


(2.23) 


where  A  Is  a  nonsingular  NxN  matrix,  Z  Is  a  Nxl  column  and  o,  0  are 


scalars  related  by 


«_1  +  0"1  •  zW 


(2.24) 


Using  the  matrix  Inversion  lemma  to  Invert  *  In  Equation  (2.22)  one 


♦'*  ■  ^  ("n1  -  *  Ud  U5  CxO*)'1 


where 


i  •  L.  ♦  uj  R-1  u! 
X  d  n  d 


(2.25) 


(2.26) 


The  array  will  acquire  and  track  the  desired  signal  If  the 
reference  signal  Is  correlated  with  the  desired  signal  and  Is 
uncorrelated  with  Interference  signals.  Assuming  that  the  reference 


signal  R(t)  Is  given  by, 

R(t)  »  Arej(°k>t  +  *d) 

and  using  Equation  (2.13),  Equation  (2.17)  yields 

,  _i.*  * 

S  *  Ap  Ad  (ZQl)  Ud  • 

Using  Equations  (2.25)  and  (2.28)  the  steady  state  weights 
(Equation  (2.15))  of  the  array  are  given  by 


(2.27) 


(2.28) 


V 


(2.30) 


K  =  d.  ( 1  —  t  U T  R'1  U*) 
~3  K  and-' 


Is  a  constant.  The  weights  given  in  Equation  (2.29)  will  lead  to  the 
maximum  output  SINR  in  the  presence  of  multiple  jammers  [see  Appendix]. 
Knowing  the  steady  state  weight  vector,  one  can  compute  the  output  SINR 
of  the  array  which  is  given  by 


SINR  = 


l  Pile  +  Pn 

k=l 


where  Pq  is  the  output  desired  signal  power. 


(2.31) 


pd  -  jEti(z;ixd)Twiz> =4  iud  oy  • 


(2.32) 


Plk  is  the  output  interference  power  due  to  the  k^  jammer 
Plk  =  i  Et|(Z-1X1k)TW|2> 


iuu  <z;‘>Y 


and  Pn  is  the  output  thermal  noise  power 


(2.33) 


Pn  =  |W|2  • 


(2.34) 
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(2.35) 


Using  Equations  (2.32)-(2.34),  Equation  (2.31)  yields 

SINR  *  Ka  uj  R"1  U*  • 
a  d  n  d 

Equation  (2.35)  Is  used  to  compute  the  steady  state  output  SINR  of  an 
adaptive  array  consisting  of  N  half-wavelength,  center-fed  dipoles.  All 
the  dipoles  are  assumed  to  have  similar  radiation  characteristics  and 
are  spaced  at  a  distance  d  apart  (Figure  3).  The  desired  signal  and  all 
jammers  are  assumed  to  be  theta  polarized  (Figure  3).  For  the  results 
presented  here,  f|<(9,4>,p)  =  1.28142  [Equation  (2.101. 

Figure  4  shows  the  output  SINR  of  an  adaptive  array  of  six  dipoles 
as  a  function  of  the  desired  signal  direction.  The  dipoles  are  spaced 
at  a  distance  of  half  a  wavelength  and  each  dipole  is  terminated  in  a 
load  impedance  equal  to  the  complex  conjugate  of  the  self  impedance  of  a 
half -wavelength,  center-fed  dipole.  The  input  signal -to-noise  ratio 
(£d)  Is  5  dB  and  the  output  SINR  is  computed  in  the  absence  of  all 
jammers.  The  continuous  curve  in  the  figure  shows  the  output  SINR  when 
the  mutual  coupling  between  the  array  elements  Is  taken  into  account 
while  the  broken  curve  ( — )  represents  the  output  SINR  when  the  mutual 
coupling  between  the  array  elements  is  ignored.  Note  that  the  presence 
of  mutual  coupling  changes  the  array  performance  and  the  output  SINR  of 
the  array  depends  on  the  angle  of  arrival  of  the  desired  signal.  The 
dependence  of  the  array  output  SINR  on  the  angle  of  arrival  of  the 
desired  signal  can  be  explained  as  follows.  Mutual  coupling  changes  the 
desired  signal  component  of  the  element  output  voltages  [Equation 
(2.6)].  The  array  illumination  due  to  the  desired  signal  is  no  longer 
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Figure  4.  Output  SINR  of  an  array  of  six  half-wavelength,  center-fed 
dipoles  vs.  the  desired  signal  direction  (<J>d).  Cd=$  d8, 
9d=90°,  d=0.5X,  ZL=Z*1  • 


uniform  and  depends  on  the  angle  of  arrival  of  the  desired  signal,  while 
the  noise  being  internal  Is  not  affected  by  the  mutual  coupling  between 
the  array  elements.  The  output  SINR  of  the  array,  therefore,  changes 
with  the  angle  of  arrival  of  the  desired  signal. 

In  the  above  example,  the  interelement  spacing  was  large  (x/2)  and 
we  found  that  the  mutual  coupling  does  nevertheless  affect  the 
performance  of  the  array.  For  small  Interelement  spaclngs,  the  mutual 
coupling  between  the  array  elements  will  be  large  (Figure  5)  and, 
therefore,  the  array  performance  will  be  affected  more.  This  is  evident 
In  the  plots  of  Figure  6,  where  the  output  SINR  of  the  array  is  plotted 
as  a  function  of  the  Interelement  spacing.  The  desired  signal  Is 
Incident  from  the  broadside  direction  (90°, 0°)  and  jammers  are  assumed 
to  be  absent.  The  broken  curve  In  the  plot  shows  the  output  SINR  in  the 
absence  of  mutual  coupling  while  the  continuous  curve  shows  the  output 
SINR  when  the  mutual  coupling  Is  taken  Into  account.  Note  that  the 
mutual  coupling  between  the  array  elements  affects  the  array  performance 
even  for  large  Interelement  spacing  (d  >  */2).  The  effect  is  more 
pronounced  for  small  Interelement  spacing  (d  <  ty2),  where  the  output 
SINR  drops  below  the  expected  value  (In  the  absence  of  mutual  coupling) 
by  a  significant  amount.  The  change  In  the  output  SINR  can  again  be 
explained  using  Equation  (2.6).  In  the  presence  of  mutual  coupling,  the 
array  Illumination  due  to  the  desired  signal  deviates  from  the  uniform 
Illumination.  The  amount  of  deviation  depends  upon  how  strongly  the 
array  elements  are  coupled.  The  noise  being  internal  to  the  processors 
Is  not  affected  by  the  mutual  coupling  between  the  array  elements.  The 
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Figure  6.  Output  SINR  of  an  array  of  six  half -wavelength,  center-fed 
dipoles  vs,  the  Interelement  spacing.  £d=5  dB,  Zl=Z-| -j * , 

( ®d»  ^d)=(90°,  0°). 


output  SINR  of  the  array,  therefore,  changes  with  the  Interelement 
spacing. 

The  drop  In  the  output  SINR  for  small  interelement  spacings  can 
also  be  related  to  the  reduction  of  the  total  incident  energy.  As  the 
Interelement  spacing  is  decreased,  the  total  aperture  of  the  antenna 
decreases  and  so  does  the  total  incident  energy  due  to  the  desired 
signal.  Since  the  receiver  internal  noise  remains  unchanged,  the 
signal -to-nolse  ratio  drops.  For  the  same  reason,  the  introduction  of 
additional  elements  into  a  fixed  aperture  with  half  wavelength  spaced 
antenna  elements  can  degrade  the  adaptive  array  performance.  The  total 
aperture  is  fixed  and  so  is  the  total  incident  energy.  The 
introduction  of  additional  elements  will  add  to  the  thermal  noise 
without  increasing  the  available  signal  power  and  that  will  degrade  the 
array  output  SINR.  Figure  7  shows  the  output  SINR  of  an  adaptive  array 
as  a  function  of  the  number  of  antenna  elements.  The  array  is  a  linear 
array  of  half -wavelength,  center-fed  dipoles.  The  total  aperture  is 
fixed  at  2 A  and  the  desired  signal  is  incident  from  the  broadside 
direction.  Again  jammers  are  absent.  The  output  SINR  is  computed  with 
and  without  mutual  coupling.  In  these  plots  only  the  indicated  points 
are  meaningful  (the  total  number  of  antenna  elements  is  always  an 
integer).  Note  that  in  the  absence  of  mutual  coupling  the  output  SINR 
increases  with  the  introduction  of  additional  elements.  It  is 
consistent  with  the  previous  work  of  Compton  [51.  But  in  the  presence 
of  mutual  coupling,  the  array  output  SINR  decreases  with  the 
introduction  of  additional  elements.  One  can  see  that  the  output  SINR 
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Figure  7.  Output  SINR  of  an  array  of  half-wavelength,  center-fed 
dipoles  of  fixed  aperture  vs.  the  number  of  elements. 
5,=5  dB,  Z,  =Z*  (9  ,,<!>.)= (90°,  0°),  total  aperture  =  2 


reaches  a  maximum  for  a  four-element  array  and  the  array  performance 
degrades  with  the  Introduction  of  additional  elements.  The  array 
performance  Improves  as  the  number  of  antenna  elements  Is  Increased  from 
two  to  three  to  four.  The  reason  for  this  Is  that  one  needs  a  minimum 
number  of  antenna  elements  to  receive  all  the  energy  Incident  on  a  given 
aperture.  But  beyond  this  point  the  aperture  becomes  overcrowded  and 
thus,  a  densely  packed  array  (d  <  ty2)  may  lead  to  a  worse  performance 
In  an  adaptive  mode. 

In  the  examples  given  so  far,  the  array  performance  was  computed  in 
the  absence  of  jammers.  The  presence  of  jammers  will  degrade  the  array 
performance.  As  pointed  out  in  our  earlier  work  [1],  the  degradation  in 
the  array  performance  can  be  computed  using  the  unperturbed  pattern  of 
the  array.  The  unperturbed  pattern  of  an  adaptive  array  was  defined  to 
be  proportional  to  the  radiation  pattern  of  the  array  responding  to  a 
single  desired  signal  in  the  absence  of  interfering  signals.  In  the 
absence  of  jammers,  the  normalized  noise  covariance  matrix  [Equation 
(2.21)]  becomes 

R  =  Z*  ZT  •  (2.36) 

n  o 


Substituting  Equation  (2.36)  in  Equation  (2.29),  the  steady  state 
vector  of  the  array  will  be 


-1 


W  =  K(Z0  Ud)*  . 


(2.37) 


Thus,  the  weight  vector  for  the  unperturbed  pattern  will  be  (Z0  Ud)* 
and  the  value  of  the  unperturbed  pattern  in  the  direction  (9,  $)  for 
polarization  p  will  be 

E(9,  ♦,  P)  =  (Z0“\j)T  (Zq"1^)*  (2.38) 

where  U  is  the  signal  vector  of  the  array  in  direction  (9,  <t>)  for 
polarization  p.  Substituting  Equation  (2.36)  into  Equation  (2.35),  the 
output  SINR  of  the  array  In  the  absence  of  all  jammers  will  be 

SINR  =  §d  (Z^Ud)1  (Zq  Ud)* 

*  5d  E(9d,  <j>d,  pd)  (2.39) 

and  is  proportional  to  the  value  of  the  unperturbed  pattern  (Equation 
2.38))  in  the  desired  signal  direction.  Further,  following  the  same 
procedure  as  given  in  reference  [11,  it  can  be  shown  that  the  output 
SINR  of  the  array  in  the  presence  of  one  jammer  will  be  given  by 

SINK  ■  -  |E(911-  ♦».  PU)f  ]  (2.40) 

dlldddJ  <z;l  Un)T(z;l  UU>* 

where  E(9-n,  p^)  is  the  value  of  the  unperturbed  pattern  in  the 

jammer  direction.  The  same  can  be  done  for  multiple  jammers.  Thus,  the 
degradation  in  the  array  performance  can  be  computed  using  the 
unperturbed  pattern. 

In  this  section,  the  effect  of  mutual  coupling  between  the  array 
elements  on  the  steady  state  performance  of  an  LMS  type  adaptive  array 


was  presented.  It  was  shown  that  though  LMS  adaptive  arrays  produce  the 
maximum  obtainable  output  SINR  (see  Appendix),  their  performance  Is 
affected  by  mutual  coupling.  One  should,  therefore,  take  mutual 
coupling  Into  account  to  compute  the  true  output  SINR  of  the  array.  In 
the  next  section,  the  transient  response  of  an  adaptive  array  In  the 
presence  of  mutual  coupling  will  be  studied. 

III.  TRANSIENT  RESPONSE  OF  AN  ADAPTIVE  ARRAY  IN  THE  PRESENCE  OF 
MUTUAL  COUPLING 

The  speed  of  response  of  an  adaptive  array  Is  controlled  by  the 
eigenvalues  of  Its  signal  covariance  matrix  [3].  As  pointed  out  In  the 
last  section,  the  presence  of  mutual  coupling  between  the  array  elements 
affects  the  Input  signals  to  the  adaptive  processor  and  thus  the 
covariance  matrix.  The  eigenvalues  of  the  covariance  matrix  will, 
therefore,  be  different  than  those  In  the  absence  of  mutual  coupling. 

In  this  section,  the  transient  response  of  an  adaptive  array  In  terms  of 
the  eigenvalues  of  Its  covariance  matrix  will  be  studied.  From 
Equation  (2.20),  the  covariance  matrix,  #,  can  be  written  as 

♦  ■  °z  ti  ♦  j,  v*  vT  +  *d<z;1  «/  <z;!  uik>TJ  • 

(3.1) 

In  the  presence  of  m+1  signals  (one  desired  signal  and  m  jammers), 

*  has,  at  least,  N-m-1  eigenvectors  (N  Is  the  total  number  of  antenna 

elements)  having  unity  eigenvalues  (assuming  and  the  rest  of  the 


eigenvectors  have  eigenvalues  larger  than  one.  The  presence  of  mutual 
coupling  between  the  array  elements  will  affect  these  eigenvalues.  We 
will,  therefore,  compute  the  nonunity  eigenvalues  to  study  the  effect  of 
mutual  coupling  on  the  transient  response  of  adaptive  arrays. 

Figure  8  shows  the  nonunity  eigenvalue  of  a  six  half -wavelength, 
center-fed  dipole  adaptive  array  in  the  presence  of  one  signal  only  (no 
jammer).  The  desired  signal  is  incident  from  the  broadside  direction 
(90°,  0°)  and  is  10  dB  stronger  than  the  thermal  noise.  The  eigenvalue 
is  plotted  as  a  function  of  the  interelement  spacing.  Note  that  the 
mutual  coupling  between  the  array  elements  affects  the  eigenvalues  even 
for  large  interelement  spacing,  but  the  effect  is  more  severe  when  the 
spacing  is  less  than  half  a  wavelength.  For  small  spacings,  the 
eigenvalue  drops  significantly  below  the  value  obtained  in  the  absence 
of  mutual  coupling  (continuous  line).  The  drop  in  the  eigenvalue 
indicates  a  reduction  in  the  speed  of  response  of  the  adaptive  array. 

In  other  words,  the  array  will  take  more  time  to  adapt  to  the  changes  in 
the  desired  signal  parameters.  The  main  feature  of  an  adaptive  array  is 
nulling  the  undesired  signals  (jammers).  The  transient  response  of  an 
adaptive  array  in  the  presence  of  jammers,  therefore,  will  be 
considered  next. 

Figures  9  and  10  show  the  nonunity  eigenvalues  in  the  presence  of 
one  and  two  jammers  respectively.  The  angles  of  arrival  of  the  two 
jammers  are  (90°, 30°)  and  (90°, -45°)  respectively  and  the  jammers  are  10 
dB  and  20  dB  stronger  than  the  desired  signal.  Again  the  eigenvalues 
are  plotted  as  a  function  of  the  interelement  spacing.  Note  that  the 
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Figure  8 


Nonunity  eigenvalue  (X)  of  an  array  of  six  half -wavelength 
center-fed  dipoles  in  the  presence  of  a  desired  signal  vs. 
the  Interelement  spacing.  Cd  a  10  dB,  (Sjj,  4><j)  =  (90°, 0°) 
★ 
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Figure  9. 


Nonunity  eigenvalues  (A)  of  an  array  of  six  half-wavelength 
center-fed  dipoles  In  the  presence  of  a  desired  signal  and 
a  jammer  vs.  the  interelement  spacing.  5d  *  10  dB, 

(«d.  ♦d)  *  (90°,  0°),  5il  =  20  dB.  (0ii,  ♦u)  =  (90°,  30°), 

* 


■11 ' 


27 


LCJGinm 


With  Mutual  coupling 
No  Mutual  coupling 


D  IN  WAVELENGTH 


Figure  10.  Nonunity  eigenvalues  (X)  of  an  array  of  six  half-wavelength 
center-fed  dipoles  In  the  presence  of  a  desired  signal  and 
two  jammers  vs.  the  Interelement  spacing.  5d  =  10  dB, 

( ed,  *d)  -  (90°.  0°),  511  =  20  dB.  512  -  30  dB,  (9u,  *u) 

(90%  30°),  (019,  ♦.,)  =  (90%  -45°) ,  Z,  =  Z*  . 


mutual  coupling  between  the  array  elements  affects  the  eigenvalues  even 
for  large  Interelement  spaclngs,  but  the  effect  Is  more  severe  for  small 
Interelement  spacing  (d  <  \/2),  For  such  spacing,  the  eigenvalues  drop 
substantially  below  the  values  that  they  would  have  had  In  the  absence 
of  mutual  coupling  (continuous  curve).  The  smaller  the  eigenvalues,  the 
longer  will  be  the  transient,  or,  the  array  will  take  more  time  to  null 
the  jammers,  which  may  be  undesirable.  The  strong  fluctuations  In  the 
eigenvalues  for  large  Interelement  spaclngs  are  due  to  the  fact  that  one 
or  more  signals  (jammer  as  well  as  the  desired  signal)  are  Incident 
from  grating  lobe  directions. 

In  this  section,  the  effect  of  mutual  coupling  on  the  transient 
response  of  an  adaptive  array  was  studied.  It  was  shown  that  the 
presence  of  mutual  coupling  between  the  array  elements  affects  the 
transient  behavior  of  the  adaptive  array  and  slows  Its  response  to  both 
desired  as  well  as  jamming  signals  for  closely  spaced  elements.  In  the 
next  section.  It  Is  shown  that  mutual  coupling  affects  the  performance 
of  Applebaum  type  adaptive  arrays  and  that  the  steering  vector  for  these 
arrays  must  be  modified  to  account  for  the  mutual  coupling  in  order  to 
maximize  the  output  SINR. 

IV.  APPLEBAUM  ARRAYS 

In  the  case  of  an  Applebaum  adaptive  array,  one  uses  a  steering 
vector  or  Initial  weights  Instead  of  a  reference  signal  as  the  control 
signal  (Figure  1).  The  steady  state  weight  vector  for  this  type  of 
adaptive  array  [6]  Is  given  by 


where  Us  Is  the  steering  vector,  G  Is  the  loop  gain  and  I  Is  a  NxN 
Identity  matrix.  Note  that  Equation  (4.1)  contains  the  signal 
covariance  matrix,  *.  As  pointed  out  earlier,  the  presence  of  mutual 
coupling  between  array  elements  will  change  *  and  thus  will  affect  the 
steady  state  performance  of  Applebaum  type  adaptive  arrays.  Eigenvalues 
of  the  signal  covariance  matrix  control  the  speed  of  response  of 
Applebaum  type  adaptive  arrays  too.  The  presence  of  mutual  coupling 
between  array  elements  will,  therefore,  affect  the  transient  response  of 
Applebaum  type  adaptive  arrays  In  the  same  fashion  as  pointed  out  In  the 
last  section. 

Next,  we  will  find  the  optimum  steering  vector  to  maximize  the 
output  SINR  of  Applebaum  type  adaptive  arrays  In  the  presence  of  mutual 
coupling.  Assuming  that  the  loop  gain  Is  large,  the  steady  state 
weights  [Equation  (4.1)]  becomes 

W  *  *-l  US  .  (4.2) 

Using  Equation  (2.25)  In  Equation  (4.2)  one  gets 

W  -  (R  _1  -  t  R  _1  U*  uj  R  _1)  Z*  U-  •  (4-3) 

^2  n  n  d  d  n  '  o  5 

If  the  steering  vector  Is  chosen  to  steer  the  beam  In  the  desired 

it 

signal  direction,  l.e.,  =  U^  then  the  steady  state  weight  vector  will 

be  given  by 


Comparing  Equations  (2.29)  and  (4.3)  one  can  see  that  the  two  are 
not  the  same,  and  thus  this  choice  of  the  steering  vector  would  not  give 
the  optimum  SINR.  If  Instead  of  the  "open  circuit"  desired  signal 
vector  (Uq),  the  complex  conjugate  of  the  desired  signal  component  of 
the  element  output  voltage  is  used  to  generate  the  steering  vector, 

1 .  e . , 

US  ’  (V1  “/  .  (4.5) 

then  from  Equation  (4.3) 

W  =  K,  lJn  R  “1  IK  •  (^.6) 

Ion  d 

where 

K1  -  V1  ‘  T  U4  Rn‘l  Ud>  '  (4-7) 

Comparing  Equations  (2.29)  and  (4.6),  one  can  see  that  the  two 
weight  vectors  differ  only  by  a  scale  factor.  The  choice  of  the 
steering  vector  as  given  In  Equation  (4.5)  will,  therefore,  lead  to  the 
optimum  performance  of  the  array. 


Figure  11  shows  the  output  SINR  of  an  adaptive  array  of  six  half¬ 
wavelength,  center-fed  dipoles  as  a  function  of  the  interelement  spacing 
for  the  two  choices  of  steering  vectors.  The  desired  signal  is  incident 
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Figure  11.  Output  SINR  of  an  array  of  six  half -wavelength, 
center-fed  dipoles  vs.  the  Interelement  spacing 

using  Applebaum  algorithm.  s  5  dB,  Z,=  , 

(0d.*d)  »  (90°,  0°). 
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from  the  broadside  direction  and  the  load  Impedance  is  equal  to  the 
complex  conjugate  of  the  self  Impedance  of  a  half -wavelength,  center-fed 
dipole.  Again  jammers  are  assumed  to  be  absent.  Note  that  when  the 
steering  vector  Is  chosen  according  to  Equation  (4.5),  the  adaptive 
array  gives  a  better  performance  for  small  interelement  distances  where 
the  mutual  coupling  Is  the  strongest.  Comparing  Figures  6  and  11,  one 
can  see  that  the  output  SI  NR  of  the  properly  excited  Applebaum  array  is 
the  same  as  that  of  the  LMS  array.  Thus,  IJ$  =  (Z0“l  Ud)*.  leads  to  the 
optimum  performance  of  an  Applebaum  type  adaptive  array. 

V.  CONCLUSIONS 

In  this  work,  the  effect  of  mutual  coupling  between  array  elements 
on  the  performance  of  adaptive  arrays  was  studied.  It  was  shown  that 
mutual  coupling  affects  the  performance  of  adaptive  arrays  even  for 
large  interelement  spaci  gs.  The  effect  Is  particularly  serious  for 
small  interelement  spacing  where  the  steady  state  output  SINR  of  the 
array  Is  significantly  lower  than  that  obtained  when  mutual  coupling  is 
Ignored  and  the  speed  of  response  of  the  array  is  reduced. 
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APPENDIX  A 


From  Equation  (2.21) 

m 

R„  -  Zo  f  ♦  E1k  'V*  U1k>*  (V1  U1k>T>  ZI  • 


Let 


Uik  •  V*  U1k 


then  Equation  (A.l)  becomes 


m 


Rn  ■  Zo  I'  +  J,  Sk  C  UikTt  zT 


*  T 

=  I  R'  z'  • 
o  n  o 


where 


Rn  *  1  *  Jj  51k  «ik*  UikT 


(A. 


(A. 


(A. 


Using  Equation  (A. 2)  in  Equation  (2.29),  the  steady  state  weight 


vector  of  the  array  is 


W  -  (Z0-‘  Ud)« 

-  K(RJ)'1  «i*  •  (A.4) 


where 


(A.5) 


Comparing  Equation  (A.4)  with  the  optimum  weight  vector  [Reference 


6,  Equation  (4.1)]  one  can  see  that  the  two  are  similar.  Thus  the 
weight  vector  given  by  Equation  (2.29)  will  lead  to  the  maximum  output 
SINR  in  the  presence  of  multiple  jammers. 
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$  MISSION  5 

1  of  I 

^  Rome  Air  Development  Center 

0  RA VC  plan*  and  executes  resear c.h,  development,  test  and  \ 

•  selected  acquisition  programs  In  support  of  Command,  Control  Q 

6  Communications  and  Intelligence  (C3I)  activities.  Technical  ? 

%  and  engineering  support  uxithln  areas  oh  technical  competence  C 
&  is  provided  to  ESV  Program  Offices  (P0<&)  and  other  ESV  * 

D  elements.  The  principal  technical  mission  areas  are  ^ 

a.  communications,  electromagnetic  guidance  and  control,  sur-  $ 

ft  velllance  oh  ground  and  aerospace  objects,  intelligence  data  r 

collection  and.  handling,  information  system  technology,  Q 

&  ionospheric  propagation,  solid  state  sciences,  mlcrom.ve 
%  physics  and  electronic  reliability,  maintainability  and 
&  compatibility. 
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